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We study holographic entanglement entropy in the background of charged dilatonic black 
holes which can be viewed as holographic duals of certain finite density states of M = 4 
super Yang-Mills. These charged black holes are distinguished in that they have vanishing 
ground state entropy. The entanglement entropy for a slab experiences a second order 
phase transition as the thickness of the slab is varied, while the entanglement entropy for 
a sphere is a smooth function of the radius. This suggests that the second scale introduced 
by the anisotropy of the slab plays an important role in driving the phase transition. In 
both cases we do not observe any logarithmic violation of the area law indicative of hidden 
Fermi surfaces. We investigate how these results are affected by the inclusion of the Gauss- 
Bonnet term in the bulk gravitational action. We also observe that such addition to the 
bulk action does not change the logarithmic violation of the area law in the backgrounds 
with hyperscaling violation. 
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1. Introduction and summary 

Understanding strongly interacting compressible states of matter has been a subject of 
intense research (see e.g. [1] for a recent review). In this pursuit, it is natural to use 
the AdS/CFT correspondence, which provides an analytic machinery of dealing with 
strongly interacting field theoretic systems. A natural route involves starting with the 
well-understood duality between the M = 4 Super Yang Mills in 3 + 1 dimensions and type 
IIB supergravity in AdS$ x S 5 and deforming it by turning on global charge density. The 
simplest setup — a truncation to Einstein-Maxwell theory — gives rise to the extremal RN- 
AdS black hole at vanishing temperature; fermion correlators in this background exhibit 
poles which can be associated with (non)Fermi liquid behavior [2-4]. 

The AdS2 near-horizon geometry of the extremal RN-AdS black hole is the bulk 
feature responsible for this exotic behavior. However, the generic extremal black hole also 
has a finite entropy density at zero temperature signalling that the groundstate is unstable. 
It is natural to ask whether there are finite density states, preserving the AdS2 factor, that 
do have vanishing entropy at zero temperature. The answer is affirmative; turning on two 
out of three diagonal U(l) global charges in the IIB AdS§ x S 5 supergravity gives rise 
to a special class of charged dilatonic black holes studied in [5-6] (see also [7] for recent 
developments). The ten-dimensional string metric is known and can be dimensionally 
reduced down to five dimensions. The resulting near-horizon geometry is conformal to 
AdS2 x i? 3 , but the conformal factor ensures that the entropy at small temperatures 
vanishes linearly, S ~ T. 

As shown in [6] , fermion correlators in this background can be computed analytically 
and exhibit the sought-for (non-)Fermi liquid behavior, but now in an entropically sta- 
ble system. In addition the system may carry additional fermionic excitations "behind 
the horizon" [8]. A potential probe of the presence of these "hidden" Fermi surfaces is 
the entanglement entropy (EE). As discussed in [9-12], in d + 1 dimensional theory en- 
tanglement entropy of a ball of radius R in the presence of the Fermi surface behaves 
like Sbaii, fs ~ -R d-1 log.R, which violates the naive area law, Sbaii ~ R d ~ x . In general, 
computing entanglement entropy in the interacting field theories is a very nontrivial ex- 
ercise. Fortunately, in the context of holographic duality, a simple prescription has been 
formulated by Ryu and Takayanagi [13-14]. Recent work on holographic matter with hy- 
perscaling violation [8] (following [15]; see also [16-25] for subsequent work on holographic 
entanglement in the presence of hyperscaling violation) identified a set of examples, labeled 
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by the dynamical critical exponent z, which exhibit such logarithmic violation. All these 
examples involve the hyperscaling violation parameter set to be 6 = d — 1 and feature the 
low temperature behavior of the thermal entropy S ~ T^ d ~ 9 ^ z ~ T x l z . 

It has been pointed out in [26] that in order for models with generic values of z and 
6 to holographically describe Fermi liquids, one should take the limit z— >oo, )> — oo 
with the limit 6/z held fixed. This is necessary in order to have a finite spectral density 
at finite values of momenta. Interestingly, the entropically stable charged dilatonic black 
holes studied in [6] precisely fall in this class. As we have already remarked, it also exhibits 
specific heat which vanishes linearly in temperature. As pointed out in [26], the simplest 
computation of entanglement of a belt in this geometry involves a phase transition and 
the resulting hypersurface in the bulk becomes disjoint when the belt becomes sufficiently 
thick. The latter phase of course signals the simple area law for the entanglement entropy. 

In this paper we compute the entanglement entropy of a ball, and observe that contrary 
to the slab case, there is no phase transition to the disjoint configuration. However, the 
numerical results indicate that in the limit of a large size, the entanglement entropy of a ball 
still exhibits the area law. The explicit absence of a phase transition for ball-entanglement- 
entropy is odd compared to the belt configuration. One may think of the belt as the limit 
of two concentric balls with large radii but fixed difference. Geometrically, the belt phase 
transition is then of the Gross-Ooguri type for concentric Wilson loops [27-32]. When the 
difference between the radii is small, the minimal embedding surface is a half-torus; when 
the difference is large, the minimal embedding is two concentric spheres. 

We also investigate entanglement entropy in the presence of higher derivative (Gauss- 
Bonnet) term in the gravity action to address how general these results are. The pre- 
scription for computing entanglement entropy in the presence of such a term has been 
formulated in [33-35], and we assume that it is not modified in the presence of additional 
matter in the bulk. In the holographic Gauss-Bonnet gravity we restrict our attention to 
the EE of a slab in the geometry with hyperscaling violation (which described the near 
horizon region the the Einstein-Hilbert case). Already here we encounter a surprise: the 
connected entangling surface does not end on the boundary for large values of the thickness 
i of the slab; it can instead approximate the boundary but never get there. We expect that 
completing the geometry would lead to the same picture as that for the ball entanglement 
entropy in the Einstein-Hilbert case: no phase transition and a connected surface in the 
bulk defined for all large values of I. 
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The paper is organized as follows. In the next Section we review both the 5- 
dimensional dilatonic black hole and its 10-dimensional lift and show that entanglement 
entropy calculation is equivalent in both settings. In Section 3 we compute entanglement 
entropy of simple geometries and observe that the calculations in the case of a belt and 
and in the case of a ball are very different. In the former case there is a maximal value of 
the belt thickness, which can support a curved entangling surface; beyond that value there 
is a phase transition to a disjoined entangling surface. In the case of a ball, the entangling 
surface in the bulk stays connected and arbitrary large values of the ball radius R are sup- 
ported. However, the entanglement entropy exhibits an area law in both cases. In Section 
4 we study EE in the case of Gauss-Bonnet gravity and comment on the differences with 
the Einstein-Hilbert case. We discuss our results in Section 5. In the appendix we present 
some details of EE in Gauss-Bonnet gravity computed in the holographic duals of states 
with hyperscaling violation. 

Note Added: When this paper was completed [36] appeared where the entanglement 
entropy of a straight belt in the five-dimensional black hole geometry was computed. This 
partially overlaps with Section 3. 

2. Entanglement entropy and dimensional reduction 

As pointed out in the introduction, one of the characteristic properties of systems with 
a Fermi surface, is the logarithmic violation of the area law of entanglement entropy [9- 
12]. Such logarithmic violations were recently found [8-16] in a class of holographic, albeit 
singular (with exception the case 9 > and z = l + ^in<i + 2 bulk spacetime dimen- 
sions 1 ), geometries with flux. It was subsequently realized that these class of geometries 
additionally exhibit other characteristic properties of Fermi surfaces, such as the scaling 
of the entropy with temperature at low temperatures, i.e., S ~ [8]. 

The geometries in question are characterized by a critical exponent z and a hyper- 
scaling violation exponent 9 and are described by 

*** = ^ (~ r 2 diz -t(d-0) + 9orT^dr 2 + £ dx 2 j , (2.1) 

1 We thank Edgar Shaghoulian for correcting a mistake in this statement in the earlier version 
of the article. 
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where the UV boundary is located at r = and the IR corresponds to large values of r. 
They can be found as solutions of the field equations emanating from 

£ = R- Z($)F 2 - (d$) 2 - V($) (2.2) 

where the dilatonic scalar field $ and the Maxwell field A are additionally turned on. The 
parameter go which appears in (2.1) is related to the IR behavior of the potential V(<&). 

This class of models has been extensively discussed in [37-41]. The logarithmic viola- 
tion of the area law in entanglement entropy is realized for 9 = d — 1. Scale transformations 
act on these geometries [8] 

it ->■ Kit 
t ->■ AH 

(2-3) 

ds -> Adds 

(d-6) 

r — > A d r 

(d-6) 

and lead to the following behavior for the thermal entropy S ~ T * , explicitly verified 
in [15] and [8]. One recognizes how the hyperscaling violation exponent theta parametrizes 
the effective dimensional reduction. 

An interesting special case of (2.1) involves the double scaling limit 

Q 

z ->■ oo = ?7>0, (2.4) 

z 

with rj fixed. As emphasized in [26] several physical quantities behave as expected for 
Fermi liquids in this limit. The spacetime metric becomes conformal to AdS2 x R d 

l9 1 ( dt 2 dr 2 J-^ , o\ 

r \ rn r i=1 J 

while the thermal entropy vanishes like S ~ T v for small temperatures. More importantly 
the spectral density remains finite at low energies and finite momenta [26] . 

Here we are interested in a specific example of (2.5), the dilatonic black hole in AdS$ 
recently explored in [5-6]. This black hole is usually referred to as the two-charge black 
hole, because it arises as a truncation of IIB supergravity on AdS 5 x S 5 where two of the 
three U(l) charges are equal and no n- vanishing while the third one is zero [42]. 

After this truncation one obtains the Lagrangian 

C = R - ±e 4a F 2 u - 12 {d^af + -L ( 8e 2a + 4e~ 4a ) ) (2.6) 
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with a a scalar field and the field strength of the Maxwell field A^. The extremal 
black hole metric 



with 



ds 2 = e 2A (-hdt 2 + dx 2 ) + 



A = ln^ + iln + % 



,2B 



h 



-dr z 



(2.7) 



B - -In-- - In [! + %■ 



h = 



L 3 

(r 2 + 2Q 2 )r 2 
(r 2 + Q 2 ) 2 



(2-f 



together with 



" L(r 2 "+Q 2 )^' 



a = - In I 1 + — 



(2.9) 



form a solution of the equations of motion coming from (2.6). In these coordinates, large 
r corresponds to the near-boundary UV region while the black hole horizon is located at 
r = 0. To be precise, r = is actually a naked singularity. Any finite temperature will 
cloak the singularity, however, so it is of the "good" type. In the near-horizon region, (2.7) 
is conformal to AdS 2 x -R 3 [6] 



ds' 



Q 



L 2 



2r2 dt 2 + ^dr 2 + ^±d X 2 

i=i 



(2.10) 



With a change of variable from r to r such that ^ = the IR metric reduces to 



2iL 2 



dt 2 9L 4 d? 2 
+ 



2Q 2 r 2 



(2.11) 



i=l 



Thus (2.11) is of the form (2.5) with g = and rj = 1. 

By the definition of consistent truncation, this solution can be lifted to a solution 
of the full IIB supergravity with Qi = Q 2 = Q and Q 3 = [42]. This ten dimensional 
solution is 



ds 2 = y/AH'^ 
L 2 H 



3. A 



-hdt 2 + dx 2 ] + 



,2B 



i=l 



h 



-dr z 



+ 



£ 

i=l 



d[i 2 + ix i dcj)i + 



+ 



(cos 2 e^O 2 + sin 2 9^1) 



r 4 s 

F 5 = G 5 + *G 5 G 5 = dB 4 B 4 = --j^Adt A d 3 x — Vid<Pi A d 3 x . 



i=l 



(2.12) 
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where 



O 2 

H = l + -\, A = if (cos 2 d\ + -ffsin 2 6>i) 



and h, A, B are given in (2.8). Here Hi and ^2 are defined as 



Hi = cos^i cos 6*2, A* 2 = cos 6>i sin #2 



(2.13) 



(2.14) 



where 6>i, 6*2 are coordinates on the S 2 inside the S^which is parametrized by (6>i, 6*2, <j>i, 4>2, 1 
The solution describes N coincident D3 branes, rotating with equal angular momentum 
in the two out of the three possible planes of rotation transverse to the D3-brane world- 
volume. 

As explained in the introduction, our main focus will be on entanglement entropy. A 
natural question to ask is whether the computation of entanglement entropy in the ten 
dimensional geometry (2.12) is equivalent with the one in five dimensions. In the simple 
case of p-branes, [43] and [17] showed that the two computations, in the ten dimensional 
setup and in the dimensionally reduced geometry, agreed as long as the former was done 
in string frame and the latter in the Einstein frame. In what follows, we will show that 
this is also true for (2.12) and (2.7). 

Let us consider the entanglement entropy of some arbitrary connected region and 
parametrize the entangling surface by a single function r(x\, £2, ^3)- The induced metric 
for a surface in the geometry of (2.12) which asymptotes to the entangling surface under 
consideration is 

D 2B-2A 1 

5ij + - — (dir)(djr) dx l dx : ' + 



ds 



2 

EE 



+ 



'AH 
L 2 



h 2A 



h 



[Hc\d9l + Hc\ (cjdcpl + s\d<$) + s\d^l] + 



(2.15) 



del 



I |_ J. A _L \ ' _L ' Zr'Z,/ 1 1 ' OJ ' 

where we introduced the notation a = cos#i, Si = sin#i and used the definition of A in 
(2.13) to simplify the expression. The entropy functional is proportional to the volume of 
the induced metric, which integrated over the five dimensional compact space yields 



y/9EE = Vol 5 e 



3A 



'det 



o2B-2A 



h 



(dir)(djr) 



(2.16) 



where V0I5 = L 5 Os is the volume of the five dimensional sphere of radius L. Finally, the 
entanglement entropy computed with the ten-dimensional metric (2.12) is given by 



qEE 
D slab 



If 



d 3 xe 3A 



'det 



8ij + dirdjr- 



D 2B-2A 



h 



(2.17) 



where we used the relation = -^nf. As expected, (2.17) precisely coincides with the 
entropy functional computed in the dimensionally reduced metric (2.7). Note that in the 
ten-dimensional metric the singularity is that of a conventional rotating black brane. The 
equivalence thus guarantees that the apparent naked singularity in the five-dimensional 
extremal black hole metric (2.7) is innocuous and does not affect the results. 

In summary, computing entanglement entropy in the ten-dimensional geometry is 
completely equivalent to computing it in the five dimensional geometry. 



3. Entanglement entropy for the dilatonic black hole 

In this Section we will compute holographic entanglement entropy in the theory described 
by the metric (2.7). In the following we set L = 1. It will be useful to write the near-horizon 
limit of the metric (2.7) in the following form 

9 5^2 2 f —dt 2 + (1Z 2 .\ . 

cfe 2 = 2-3Q-^-3 ^ _ + J2 dx l) C 3 - 1 ) 

where x^ = y/^Qx^ and z = l/2r. Note that the metric (3.1) belongs to the class of metrics 
studied in [8] with a particular value of dynamical exponent z = oo. We start by placing a 
system in the large 3-dimensional box of size L x , which serves as an IR regulator. We can 
now compute the entanglement entropy of a slab, specified by < x\ < £, < x%, xs < L x , 
where I <C L x . The entangling surface can be parameterized by a single function x\(z) 
with the boundary conditions xi(0) = 0, 1. The entropy functional takes the form 



S slab = L 2 x J ^l + (x') 2 (3.2) 

where x = \p2Qx\, the prime denotes derivative with respect to z and the integral needs 
to be appropriately regularized. In (3.2) and in the rest of the paper we omit an overall 
numerical factor, together with a factor of L 2 /£ 2 , which roughly measures the number of 
degrees of freedom (in the case of the holographic dual of a conformal field theory such 
a factor is proportional to the central charge). The equation of motion that results from 
(3.2) can be written as 

T> 1 

(3.3) 



Zy/Z~ 2 + (X') 2 Z 
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where 1/zq is a constant of motion; the value of z = Zq corresponds to the turning point of 
the curve in the bulk where x' = oo. By varying the value of z$ one can change the shape 
of the curve. The value of £ as a function of z can be computed via 

/C z ° dz 1 
dzx' = 2 / = 7T. (3.4) 

Jo Z y/(zo/z) 2 - 1 

The key feature of this result is its independence of the value of z$: the curved hypersurface 
in the bulk which approximates to the boundaries of the slab on the AdS boundary can 
only do so for a fixed thickness of the slab I. We are, of course, interested in the behavior 
of the EE as £ is taken to be large. At this point it is worth recalling that there is always a 
trivial solution x' = 0: this becomes the only available solution for Q£ ^ tx/\/2. It is clear 
that the entanglement entropy for such a solution is independent of £. (this feature of the 
geometry (2.7) has been pointed out before; see e.g. [26].) 

When the full geometry (2.7) is considered, we might expect the picture to be the 
following. For sufficiently large £ the hypersurface in the bulk should probe the infrared 
limit of the geometry (2.7), described by (3.1). Therefore the discussion of the previous 
paragraph applies: there is a maximal value of £ = £ cr u which corresponds to the curved 
solution; beyond it only the trivial x' = solution exists. The transition between the 
two configurations is second order, since the curved solution asymptotically approaches 
the trivial one as £^r£ C rit- For sufficiently small £ the entangling surface probes the UV 
(small z region) of the geometry, and we expect £ to be a smooth function of zq with a 
limit £(z = 0) = 0. To check these expectations we repeat the calculation above for the 
full geometry (2.7). The EE functional is now given by 

S slab = L 2 x j dze 2A y/e™/h + e* A (x>y (3.5) 

The analog of (3.4) is now given by 

r°° dr 
Q£ = 2 (3.6) 

J -o yV 2 (2 + r 2 ) (r 2 (l + r 2 ) 2 r - 2 (l + r 2 )" 2 - l) 

where r = r/Q, r = r /Q. The value of this integral as a function of r is plotted in Fig. 
1. 
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Fig. 1: IQ as a function of tq computed using eq. (3.6). 




Fig. 2: AS as a function of IQ. 

The figure is in accord with the general discussion above. In the UV region, tq ^> 1, 
the length of the interval monotonically decreases to zero. On the other hand, there is a 
limiting value of £ as the position of the tip approaches the horizon at r = 0. Note that 
this limiting value is equal to £ cr u = 7r/y2Q, in accord with the near-horizon analysis. We 
check that for £ < £ cr a the curved solution is preferred (gives the leading contribution to 
the EE), and the difference vanishes in the limit £^-£ C rU- (see Fig. 2). 

The next step is computing the entanglement entropy for a sphere, defined by 
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We start by considering the near-horizon metric (3.1). The entangling surface in the bulk 
can be parameterized by a single function z(p); the EE functional is given by 

S spher e = J %V(^)7^ + 1- (3-8) 

The equation of motion is now second order: 







p 2 z' 



p 2 (z 2 + 2(z') 2 ) 



(3.9) 



dp \z 3 y/(z') 2 /z 2 + l J Z 4 yS(z') 2 /z 2 + 1 

We can integrate this numerically; the result is given by the red curve in Fig. 3. The 
surprising feature of eq. (3.9) is the absence of the solutions approaching the boundary at 
z = at finite p. Such a behavior has not been observed in the examples studied in [8]. 
This behavior is evident in Fig. 3. (the red curve never intersects the z = line). Also, 
unlike the slab case, there are no trivial solutions p(z) = const. 

The near-horizon approximation is clearly not reliable and we must we analyze the 
EE of a sphere in the full diatonic black hole metric (2.7). The EE functional now takes 
the form 



S sphere = f %vWe^) + 



e 2B{z)( z ,)2 



*m (3 ' 10) 

where A(z),B(z),h(z) are determined by (2.8) with the substitution z = l/2r. The 
equation of motion for z(p) that follows from (3.10) is second order and has to be dealt 
with numerically. A representative solution is shown in Fig 3 (blue curve). The red curve 
corresponds to the solution in the near horizon geometry. The modification of the EOM 
in the full geometry ensures that the hypersurface now reaches the boundary at finite p. 




Fig. 3: z(p) for the near horizon (red) and full BH metric (blue) 
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/ 10 15 20 25 30 35 40 

Fig. 4: R as a function of zq for the full BH metric 

We can now compute the radius of the entangling surface at the boundary as a function 
of zq numerically. The result is shown on Fig 4. Finally we are in position to compute the 
value of entanglement entropy. (Similar computations have recently appeared in [44] and 
[45] ; we will follow the procedure outlined in [32] , which the reader is encouraged to consult 
for technical details.) The entanglement entropy of a ball is UV divergent and a suitable 
regularization is necessary. The UV divergent terms are governed by the asymptotically 
AdSs metric; these terms in the conformal case have been computed in [14] : 

R 2 1 

Sball,divergent = + ^ e (3-H) 

here e 1 is the cutoff on the radial coordinate z. An extra factor of 1/2 in the first term 
in (3.11) compared to eq. (7.10) in [14] appears due to the factor of two in the definition 
of z in our paper. 

All we need is to evaluate (3.10) on the solution of equation of motion, and subtract 
the UV divergent term (3.11). We are interested in the deviation of the resulting finite 
quantity, S uv finite = S ba ii - S baUj divergent from the area law. In Fig. 5 we plot the 
normalized UV finite part of the entanglement entropy. One observes directly that at 
large R the behavior of entanglement entropy of a ball is governed by the area law. The 
equation of motion does not appear to have a second distinct class of solutions and the 
result implies that the ball-entanglement entropy does not experience a phase transition 
at smaller radius to a phase differing from the naive area law. 

11 



0.20 - 




-0.10 - 



Fig. 5: Suv finite/ R as a function of RQ; red curve corresponds to Suv finite /Q R = 
0.25 - 0.78/Q 2 R 2 fit. 

4. Gauss-Bonnet and the EE of a belt 

Here we compute the entanglement entropy of a belt in Gauss-Bonnet gravity. We use the 
near horizon geometry of eqs. (3.1) which is a solution to the equations of motion following 
from dilatonic Gauss Bonnet gravity with flux (see the appendix for more details). En- 
tanglement entropy in Gauss-Bonnet gravity is given by extremizing the action functional 
[34-35] 

S = J d^VdetGs (1 + Alfe) (4.1) 

where A is the Gauss-Bonnet coupling constant whereas Gy, and Ry, refer to the induced 
metric and the induced scalar curvature 2 . 

Consider the slab, specified by < x\ < I and X2, £3 < L x with L x » I as in section 
3. The entangling surface is parametrized by a single function x\{z) and gives rise to the 
induced geometry 

( —z + xJz)) dz 2 + \^ dxj 

^ / u 

where the dot implies differentiation with respect to the variable z and Xi = \plQxi. The 
action functional (4.1) reduces in this case to 

dz I 1 +i 2 (l 2XZ ^ i 1 + 7Z<2±2 + 6z3± x) 



ds =2 3Q 



(4.2) 



Sslab — 2 S Q 2 L 2 



9(1 + z 2 x 



2~2\2 



(4.3) 



2 To make the variational problem well-defined a boundary term should be added to (4.1). 
This term only affects the leading (proportional to A) divergent term in the entanglement entropy, 
e.g. see [35], it will therefore not affect the results of this section. 



12 



where we set x = \[2Qx\. The equation of motion folowing from (4.3) takes the form 

x 2 2 x 1 

-n A ^ 3 : r-^r = - ( 4 - 4 ) 



VI + 2 2 x 2 9 (1 + ^2)2 z 
where z$ is a constant of motion related to the width of the slab £ = -J^q through 

I f° 

- = I dzx (4.5) 



«0 



It is convenient to define a dimensionless variable y = which takes values in y e [0, 1] 
and express (4.4) as follows 

>|- ^—^-1 = 0, (4.6) 



where now primes indicate differentiation with respect to the new variable y. Eq. (4.6) 
shows that the presence of the Gauss-Bonnet coupling strongly affects the physics of the 

2 

entanglement entropy. In particular, x'(y) now depends on zq through B = Xzq . As a 
result the length of the slab £, which in y coordinate is given by 

£ - = J\yx'{y), (4.7) 

is not anymore independent of zq. Note that we have taken zq > to make contact with 
the A = limit. 

To proceed we need to solve (4.6) for x'(y). Obvisouly, the disconnected solution, 
x' = is always a solution. There also exist three connected solutions to (4.6) but, as 
commonly observed in Gauss-Bonnet gravity, only one of them asymptotes to the connected 
solution when A vanishes 3 . We choose this solution for x' and proceed to examine its 
behavior close to the boundary y = 0. We find that for all B > | x'(y) diverges as y 
approaches the boundary 

2y -27 + 6B 2 



x'=< 



2 



y 3 A o Z 



It would be interesting to study all possible solutions in detail in the spirit e.g. of [46]. 
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whereas it goes to a constant for B < | 



x' = — — + -3* 3 + 2B 4 y 2 + 0(y 4 ) B < - (4.9) 
9-2B 2 (_ 9 + 2S) 4 



Hence, a solution which asymptotes to the entangling surface only exists for B < |. A 
natural question to ask is whether this result changes when one considers the full solution 
and not just the near horizon geometry. We leave this to future investigation. 

It is instructive to also consider the behavior of the solution close to the turning point 
y = 1. One essentially finds that the solution becomes complex for B < — |, in other 
words it ceases to exist. For B > —9/4 x' diverges at y ~ 1 and the solution tends to a 
constant value. We will therefore restrict to values of B lying in the rgion — ~ < B < |. 

Substituting this solution into (4.7) and performing the integration numerically shows 
that for B E (— |, |), the length of the slab £ varies from £ E (0.4, 3tv). However, since 
zq > 0, the sign of B is completely determined by the sign of the Gauss-Bonnet coupling 
A which according to the null energy condition (see appendix) is bounded above by A < |. 
This implies that both signs are allowed for B. It would be interesting to see if other 
constraints, such as those discussed in [47-50], can uniquely fix the sign of A. 

The length £ as a function of B for positive values of B is plotted in Fig. 6. 

4.5 - 

4.0 - 

3.5 - 

3.0 - 

2.5 - 



Fig. 6: I as a function of B = Az 3 for B > 0. A solution only exists for £ < 3n. 

Fig. 7 instead shows £ for negative values of B. The two curves join at the point B = 
where the result of section 3 is recovered; the connected solution exists only for £ = n. 
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Fig. 7: I as a function of B = Xz§ for A < 0. 

In summary, for negative A the connected solution exists for I G [0.4, tt) whereas for non- 
negative A, it exists for I G [tt, 3tt). 

We conclude this section by computing the difference between the entanglement en- 
tropy of the connected solution and the disconnected one. The difference is given by 4 
32Q 2 _32Q 2 

Y2 L ^ J bclt — T~2 \ belt, con. °belt, disc.) — 

A \ * f / . 2B\ 



d ( ( 9 ~ 2B V ls ) 2 ~ 54 (- 3 + 2B3/I ) y 2 x' 2 - 9 (-9 + 8Byf ) yV 4 _ 9 + 2 By 2 ^ 



x 9y Vl + y 2 x' 2 (9 - 2Syl + (9 + 4By§) ?/V 2 ) V y 

(4.10) 

where we used (4.6) and the fact the (4.3) is independent of x, to express the integral in 
terms of x' alone. In Fig. 8 we plot AS for positive B and in Fig. 9 for negative B. 
It is easy to see that for positive Gauss-Bonnet coupling, the connected solution-whenever 
available - is favoured, while the contrary happens for negative values of the coupling. 



5. Discussions and puzzles 

In this paper we computed holographic entanglement entropy in the background which 
corresponds to the finite density state in the Af = 4 SYM theory with two equal U(l) 

4 In the following we express the difference in the entropy in terms of B instead of t for reasons 
of convenience. Note that B is a monotonic function of I. 
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32 Q 1 AS 




Fig. 8: AS as a function of B for B > 0. 



32 g 2 AS 
Lr 




Fig. 9: AS as a function of B for £ < 0. 

charges. The computations were done for two simple geometries: a belt and a sphere. 
The puzzle we encounter is that we observed a phase transition in the case of a belt as a 
function of its width, whereas the sphere does not. 

Such phase transitions in entanglement entropies have been observed before in con- 
fining geometries [51] [43] [32]. (Of course, to make intelligent statements about the phase 
transitions, one needs to define a UV-finite quantity; one suitable definition in our case 
would involve subtracting the entanglement entropy in the vacuum, zero density state.) In 
[43] a field theoretic explanation of this effect was proposed: the value of entanglement for 
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the free glueball of mass m is exponentially suppressed, S ~ exp(— m£); integrating this 
over the glue ball density of states with Hagedorn behavior, p(m) ~ exp(/3ffw) produces 
a phase transition at £ ~ (3h- The transition, which resembles Hagedorn phase transition, 
is between the 0(1) values of entanglement for sufficiently large values of £ and 0(N 2 ) 
values for small values of £, where the integral over glue ball states formally diverges. 

One would argue that the phase transition is a universal phenomenon only depending 
on the size of the entangling surface compared to the density. Indeed in [32] a similar 
phase transition has been observed for the entanglement entropy of a sphere (disk) in 
confining geometries dual to the 3+1 (2+1) dimensional theories. This raises a second 
puzzle, since the (UV-finite part of) the free field entanglement entropy in this case is not 
exponentially suppressed, but rather has an expansion in the inverse powers of the radius 
(see e.g. [52-53]). 

We do not have a resolution of this second puzzle. Nor do we have a fundamental 
answer to the first. Instead, we add the observation that an understanding of these phase 
transitions must be more complicated in that a third scale must play a role. This third 
scale is supplied by the anisotropy of the belt. This must be so as one ought to be able 
to obtain the belt shape by deforming the ball-entangling surface continuously into an 
ovoid. In this process the phase-transition should suddenly appear; qualitatively it must 
therefore depend on the eccentricity of the ovoid. To check this explicitly is technically 
quite involved. One can, however, easily construct a different geometric set-up which 
interpolates between the ball and the belt and thereby argue that this interpretation is 
correct. This set-up is the annulus or two concentric balls. In the limit where the inner 
radius vanishes one has the sphere, and in the double limit where both radii are large, 
but the difference stays small, one has the belt. It is well-known from the analogue study 
of holographic Wilson loops that this annulus system has a phase transition between a 
geometric embedding surface that is a half-torus for the belt-like configuration and two 
nested sphere-like embeddings for the ball-like configuration [27-31]. This suggests that 
the belt-like phase transition may be a consequence of the third scale introduced by the 
anisotropy violation. It would be interesting to construct a free-field model which would 
reproduce this behavior, although there is no guarantee such a free field description exists. 

Note, that for the case of states with Fermi surfaces, one expects the two geometries 
to produce the same violation of the area law (and this is also the case in the holographic 
models of [8]). Of course, the area law is not violated in our case, so one can conclude that 
the ground state dual to the geometry explored here is different from a Fermi liquid. 
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We would finally like to note that the addition of the Gauss-Bonnet term to the 
bulk Lagrangian modifies the picture. In this case there are strong indications that the 
belt geometry can support connected surfaces for arbitrarily large values of the distance 
between the plates. Generally, adding higher derivative terms such as the Gauss-Bonnet 
term amounts to the inclusion of l/gy M N and/or 1/N corrections in the dual field theory. 
It suggests then that the phase-transition is a strong coupling/large N phenomenon which 
would be hard to capture in any field theoretic approach. The Gauss-Bonnet term is very 
special and we should be careful to extrapolate generic lessons from it. 
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Appendix A. Gauss-Bonnet, hyperscaling violation and entanglement entropy. 

A.l. Solutions of GB gravity +matter which give rise to holographic hyperscaling violation. 

Here we show that geometries with hyperscaling violation are solutions of Gauss-Bonnet 
gravity with matter and gauge fields. We consider the following five-dimensional action 

12 AL 2 
L 2 + 



(A.l) 



where A is the dimensionless Gauss-Bonnet coupling and C m is the matter Lagrangian 

CT = Z($)F 2 + -(<9$) 2 + y($). (A.2) 

Here C(F 2 ) is an arbitrary function of the quadratic gauge invariant F 2 = F^ V F^ V , with 
F^v = d^A v — d^A^ the electromagnetic tensor field and the vector potential. $ is 
a scalar field and V(<&), Z($) are arbitrary functions. The action Iq should in principle 
be supplemented with a boundary term to make the variational problem well-defined. For 
our purposes the detailed structure of the boundary term will not be necessary. 
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Varying the action with respect to the metric tensor and the matter fields <£>, 
yields the following equations of motion 

G -—a +—G {2) -T m 

V„Z(«) (F"") = (A.3) 

where = — tj-R<7 m is the Einstein tensor, T™ is the energy momentum tensor for 
the matter fields and G$ is the second order Lovelock tensor 

1,2 2 2X ^ 

The matter energy momentum tensor has the form 

= -\g»X m + 2(d^)(d l/ <S>) + 2Z{®)F£F vp . (A.5) 

We search for solutions of the form 

ds 2 = e 2A ^ (dr 2 - e 2B M >dt 2 + dxf) 
A t = A t {r) Ai = 0, i = l,2,3, $ = $(r) . 

Substituting the ansatz (A. 6) into the equations of motion (A.3) yields 



(A.6) 



T t, m _ T i,m = 2F tr F tr Z($) = -e~ 2Ai - r) (3A'(r)B'(r) + B'(r) 2 + B"(r)) + 

+ 2AL 2 e- 4A(r) A'(r) [B'{r) (A'(r) 2 + A'{r)B'{r) + 2A"{r)) + A'{r)B"{r)] 

T r,m _ T t,m = e _2A(r) = 

= 3e~ 2A ^ (l - 2\L 2 e- 2A ^A'{r) 2 ^ (A'(r) 2 + A'(r)B'(r) - A"(r)) 
8 r (e 5A+B Z($)F tr ) = 

e - 5A{r) -B(r) d ( e 3A(r)+B(r)Q Q\ = ^_ + 2 — F tr F tr . 

(A.7) 

We can easily solve the Gauss Law constraint for Z(Q)F tr in the third line of (A.7) to get 

Z($)F tr = Qe~ 5A{r) - B{r) (A.8) 
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where Q is the total charge of the solution. Substituting the result into the first line of 
(A. 7) determines F tr in terms of A(r), B(r) and their first and second derivatives. To be 
specific 

3A(r)+B(r) 

Ftr = ^ (3A'(r)B'(r) + B'(r) 2 + B"(r)) + 

^ (A 9) 

p A{r)+B{r) 1 ^ > 

+ 2AL 2 - A' [B' [A 12 + A'B' + 2A") + A'{r)B"{r)\ 

We can then use the solution (A. 9) to solve for Z($(r)) from (A. 8), 

2Q 2 e" 4A M 



(3A'B' + B' 2 + B") - 2\L 2 e- 2A ^)A'(r) [B'(r) (A'(r) 2 + A'(r)B'(r) + 2A"(r)) + A'B' 

(A.10) 



In addition, note that the second line of (A. 7) gives a solution for <9 r $ in terms of the 
metric functions A(r)B(r), i.e., 

$'(r) = ^3 (1 - 2XL 2 e- 2A ( r )A'(r) 2 ) (A'{r) 2 + A'{r)B'(r) - A"{r)) . (A.ll) 

Eq. (A.ll) can in principle be integrated to give an expression for $(r) in terms of the 
radial variable. Then it can be inverted to express r in terms of $ and determine Z(<E>). 
Finally, substituting (A.ll) together with (A. 9) and (A. 10) into the last line of (A. 7) an 
expression for ^ in terms of A{r), B(r) can be easily determined. Simply rewrite the 
derivatives in terms of $ in (A. 7) as ^ = $7^y^: and solve for ^ in (A. 7). 

It is clearly not very difficult to produce solutions of the type (A. 6) from the action 
(A.l). Not all of these solutions however are necessarily consistent either as gravitational 
solutions or as dual descitpions of quantum field theories. Absence of singularities, appro- 
priate fall-off conditions for the fields, etc. restrict the set of consistent solutions. One of 
the simplest consistency checks, is the null energy condition 

T>V>0, (A.12) 

where is an arbitrary null vector of the spacetime (A. 6), i.e., g^u^u" = 0. For the 
metric in (A. 6) the null energy condition essentially amounts to requiring the first and 
second line of (A. 7) to be non-negative 

T/' m -T/' m <0 T r r ' m -T/' m > 0. (A.13) 



20 



We will now focus on a specific choice of metric which belongs to the class of metrics 
described by (A. 6) and violates hyperscaling (for more details on these solutions and their 
condensed matter applications see [8-21]) 

ds 2 = r- 2 ^/ 3 (dr 2 - r~ 2 ^dt 2 + dx 2 ) . (A.14) 

The starting point will be to analyze the null energy condition (A. 13). Substituting the 
expressions for A{r), B(r) from (A.14) in (A. 13) yields 



(9 - 3) =2e f 2 



-26 



r~ (3 + 9 -3z) (I- -\L 2 (6 -3) 2 r— )>0 

, -26 ( 2 ~ -28 (9 -3) 2 (9 -6-3z)\ 
{z _ m _ 9 + z)r — h + -XL 2 r— y - + z ' " ° • 



(A.15) 



To determine the allowed range of values for the parameters (6, z) from (A.15) we make 
the following assumptions: 

1 that (A.14) is meaningful for all r, ranging from zero to very large values and 

2 that the theory is also consistent at the A = point. 
Eq. (A.14) then reduces to 

(0-3)(3-3z + 0) > A<0 

w n x 6-6-3z (A-16) 

(z-l)(3-e + z)>0 ____<o 

The left column of (A. 16) provides the constraints on 9, z obtained form the null energy 
condition in the case of Einstein-Hilbert gravity. It is easy to see that the right column does 
not yield additional constraints on 6, z. It only restricts the value of the Gauus-Bonnet A 
coupling to be negative. 

It is clear that the geometry of (2.5) produced by taking the double-scaling limit (2.4) 
remains a solution of Gauss-Bonnet gravity. The null energy condition in this case yields 
XL 2 < §. 

A. 2. Holographic entanglement entropy and hyperscaling violation in GB gravity 

In the following we compute holographic entanglement entropy in the GB gravity. We 
restrict to 3 + 1 dimensional boundary (so that the GB term does not vanish in the bulk). 
We will consider the general form of the metric (A. 6). The example of a metric that 
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violates hyper scaling is given by (A. 14). The holographic entanglement entropy is equal 
to the extremal value of the action 

S = J c^aVdetCs (1 + XR S ) (A.17) 

where Gy, is the induced metric. 

We consider the case of a slab as defined in section 2, with induced metric 

dsl = r- 2 ( 3 -*V 3 ((1 + (x') 2 )dr 2 + dx\ + dx 2 ) (A.18) 

where we defined x = x\. In this case the action (A.17) reduces to 

s = viTTFF) (i- 2Ac " 4,r>|j4 ' (r>2(1 + ^±^m^h^mm 

(A. 19) 

One can solve for x' and expand near the boundary x = 

-' 2 - VV + " (A ' 20) 

For < ^ < 3, x' 2 does not contribute to the UV (small r) divergence of (A. 19). 

The most interesting divergence in (A. 19) comes from the e 3A ^ term and gives rise 
to the S ~ (e/£) 2 behavior, similarly to the Einstein-Hilbert case. The case 9 = d — 1 = 
2 gives rise to a logarithmic term, S ~ log£/e, just as it did in the Einstein-Hilbert 
(A = 0) case. Interestingly, in the Gauss-Bonnet case there are terms proportional to A 
which exhibit stronger divergence, S ~ (e/£)^~ 2 . Such terms, however, do not produce 
logarithmic violations of the entanglement entropy. Moreover, they are expected to vanish 
when the contribution from the boundary term neglected in (A.17) is included. 
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